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Abstract 

This investigation extends earlier studies of a shear-transformation-zone (STZ) theory of plastic 
deformation in amorphous solids. My main purpose here is to explore the possibility that the 
configurational degrees of freedom of such systems fall out of thermodynamic equilibrium with the 
heat bath during persistent mechanical deformation, and that the resulting state of configurational 
disorder may be characterized by an effective temperature. The further assumption that the 
population of STZ's equilibrates with the effective temperature allows the theory to be compared 
directly with experimentally measured properties of metallic glasses, including their calorimetric 
behavior. The coupling between the effective temperature and mechanical deformation suggests 
an explanation of shear-banding instabilities. 
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I. INTRODUCTION 



This is the third in a recent series of studies of shear-transformation-zone (STZ) models 
of plastic deformation in amorphous solids. In the earlier two papers, Falk, Pechenik, and 
I showed how to use principles of symmetry and energy balance to constrain the form 
of STZ theories at low temperatures [if, and then used those ideas to construct a finite- 
temperature theory |2| whose predictions could be compared with the behavior of metallic 
glasses observed by Kato et al.j.3] and Lu et al. |4J . Our version of STZ theory was introduced 
originally in SL It differs from the flow-defect theories of Cohen, Turnbull, Spaepen, Argon 
and others 6,1?], Is,^ primarily in that, instead of simply postulating an equation of motion 
for the STZ density, we included a rudimentary model for the irreversible, internal dynamics 
of these zones. This augmented STZ theory exhibits an exchange of dynamic stability 
between jammed and flowing states at a stress that we identify as a yield stress. Our main 
conclusion in [al was that the transition between linear Newtonian viscosity and nonlinear 

nn 

superplasticity reported in (3, |4|] can be explained quantitatively as a transition between 
thermally activated creep at low stress and the onset of plastic flow at the STZ yield stress. 

My first purpose here is to address several issues that were left outstanding in Q. The 
most important of these were questions pertaining to the STZ density. In order to retain 
an essential simplicity in Q], we found it best to leave the theory in a form in which that 
density approached a temperature-independent value in the limit of small but nonvanishing 
driving force. We pointed out that this limiting density, after indefinitely long aging, ought 
to relax to its thermal equilibrium value, and that we would need to incorporate such an 
aging mechanism into a next version of the theory. More generally, we argued that the subtle 
interplay between limiting behaviors at small strain rates and small temperatures provides 
an important clue about the fundamental nonequilibrium properties of these systems. 

My present hypothesis is that the STZ density is governed by an effective temperature, 
T P ff, of the kind discussed, for example, in papers by Ono et al. [l^, Cugliandolo et al. 



111. Sollich et al. 1121, Berthier and Barrat 1131 . and Lacks 



0. 



Some aspects of these ideas 



are related to work by Mehta and Edwards [15|. As proposed in Refs. [1 
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T e ff differs from the ordinary thermal temperature T in circumstances where the slowly 
changing configurational degrees of freedom of the system fall out of equilibrium with the 
heat bath - a situation that occurs when molecular rearrangements are driven by plastic 
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deformation. My ideas have emerged in part from discussions with Anael Lemaitre ^(|, who 
has approached the concept of effective temperature from a different point of view. 

More specifically, I assume that, in a nonequilibrium system, the STZ density is driven 
toward n M exp (— 1/%), where x — ^BT e ff/E z . In this regard, the reduced effective tem- 
perature x is ver Y nearly, but not quite, the same as Spaepen's reduced free volume. [?], Q 
Here, Ez is a characteristic energy associated with STZ formation, and density of 

the order of the number of molecules per unit volume. This assumption implies that the 
local energy (or density) fluctuations of the slowly varying configurational degrees of free- 
dom are described by a Boltzmann distribution with effective temperature T e ff. That is, I 
assume that persistent deformation accompanied by molecular rearrangements produces a 
steady state of disorder in an amorphous system. In the absence of constraints other than 
number and energy conservation, the statistical distribution of density fluctuations must 
maximize an entropy, and therefore that distribution should be described by a temperature. 
To be consistent with the underlying assumption that the STZ's are sparsely distributed 
in the material, they must account for only a small fraction of the configurational degrees 
of freedom. Thus the probability of finding STZ's should be accurately proportional to a 
Boltzmann factor, exp (— Ez/ks T e ff), and should be far out in the wings of this statistical 
distribution. 

The underlying structure of an STZ theory based on the effective-temperature hypothesis 
is outlined in Section II. Here I reintroduce the fully nonlinear STZ theory pj in a version 
suitable for use at nonzero temperatures. As we pointed out in the so-called "quasilinear 
theory" used in the preceding papers has serious shortcomings, specifically, lack of realistic 
memory effects and unrealistically large plastic deformation at small stresses. It also has 
an unattractive asymmetry between the rates of shear transformations and dilational rear- 
rangements, which must be corrected in order to deal systematically with thermal effects. 
Section II concludes with a statement of the fully nonlinear equations of motion for the inter- 
nal state variables. In Section III, I present arguments in favor of the effective-temperature 
hypothesis and make a rough estimate of the limiting value of x a t small but nonvanishing 
strain rates. I then show how these arguments can be used to determine an equation of 
motion for x- 

The remainder of this paper is devoted to exploring the predictions of these equations of 
motion and comparing them with the results of metallic glass experiments, primarily those of 



3 



Lu et al. [4| on bulk amorphous Zr 4 i. 2 Tii 3 . 8 CU12.5 Nijo Be 2 2.5- I discuss a wide range of such 
measurements including the steady-state stress versus strain-rate data and the stress-strain 
curves obtained at various strain rates and temperatures, in analogy to our presentation in 
P). With this version of the theory, I can go on to compute specific heat curves obtained by 
differential scanning calorimetry, and also can discuss the way in which those measurements 
may be interpreted in terms of effective temperatures. I conclude with some remarks about 
how shear banding instabilities may arise in theories of the kind introduced here. 



II. BASIC EQUATIONS OF MOTION 

Let us start by summarizing briefly the assumptions and definitions used in IlUa]. Assume 
that, instead of being structureless objects as in the flow-defect theories of (J Q, , 0] , the 
STZ's are two-state systems which, in the presence of a shear stress, can transform back 
and forth between just two different orientations. Importantly, these STZ's are created and 
annihilated during irreversible deformations of the material. As in Q], consider first a two- 
dimensional system, and subject it only to pure shear deformations. (The transformation of 
the two-dimensional results into a form suitable for analysis of three-dimensional experiments 
is discussed at the beginning of Section IV. It is described in more detail in Q.) In this 
case, we need to consider only situations in which the orientation of the principal axes of 
the stress and strain tensors remains fixed. That is, we do not need to consider situations in 
which a fully off-diagonal tensorial version of the STZ theory is necessary, as in the necking 



calculations reported in 17(. Therefore, it is sufficient to assume that the population of 
STZ's consists simply of zones oriented along the two relevant principal axes of the stress 
tensor and, without loss of generality, to let the deviatoric stress be diagonal along the 
x, y axes. Specifically, let s xx = —s yy = s and s xy = 0. Then choose the "+" zones to be 
oriented (elongated) along the x axis, and the "— " zones along the y axis; and denote the 
population density of zones oriented in the "+/— " directions by the symbol n±. 
With these conventions, the plastic strain rate is: 

& = = & = ^ n - - n+) • (2.1) 

Here, A is a material-specific parameter with the dimensions of volume (or area in strictly 
two-dimensional models), which must have roughly the same order of magnitude as the 
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volume of an STZ, that is, a few cubic or square atomic spacings. The remaining factor on 
the right-hand side of Eq. (j2.1j) is the net rate per unit volume at which STZ's transform 
from "-" to "+" orientations. R(s)/t and R(-s)/r are the rates for "+" to "-" and "-" 
to "+" transitions respectively. The dimensionless stress is s — s/p,, where /2 is an effective 
shear modulus that will turn out to be an accurate approximation for the yield stress at the 
temperatures of interest here. To sets a time scale for the molecular rearrangements. As we 

n 

shall see, To is not defined here in quite the same way as it was in [2J. 

A basic assumption in this paper is that, in contrast to Eq.(3.3) in we can rewrite 
the master equation for the densities n± in the form: 

r ji± = R{t~s) n T - R(±s) n± + (T + p) e' 1 '* - n±) . (2.2) 

The first pair of terms on the right-hand side describes the same switching back and forth 
of the STZ's that appears in Eq. (j2.1|) . and the last terms describe the rates of creation and 
annihilation of zones. In writing the latter terms, I am using the principle of detailed balance 
to fix the ratio of the annihilation and creation factors, and accordingly am omitting the 
quadratic term that we used in [21]- As before, the rate factor multiplying the annihilation 
and creation terms consists of the driven part T and the spontaneous thermal part p{T). 
Our usual notation is: 

A= n+ + H -, A = n+ - U -; (2.3) 

and 

S(s) = ~(R(-s)-R(+s)), C(s) = l(R(-s)+R(+s)) : T(S) = ^. (2.4) 
Then, using Eq. ()2.1|) . and defining e = Xn^, we have 

T e pl = e C(§)(AT(§)-A); (2.5) 
r A = 2 C(s) (A T(s) - A) - (T + p) A; (2.6) 



and 

r A = (r + p)( e - 1 ^-A). (2.7) 

The next step is to use the energy-balance argument introduced in |1] to evaluate the 
quantity V. Both the effective temperature and the fully nonlinear rate factor R{s) will 
introduce features that were not present in 

Q or 0; therefore it will be useful to rewrite 



this analysis. As in the earlier papers, we start by writing the first law of thermodynamics 
in the form: 



2e pl s = C(s) (A T(s) - A) s = 4- * (A, A) + Q(s, A, A). 
r v dt 



(2.8) 



The left-hand side of Eq. (|2.8j) is the rate at which plastic work is done by the applied stress 
s = ps. On the right-hand side, \1/ is the recoverable, state-dependent, internal energy 
associated with the STZ's. Because the STZ's in this formulation represent only a very 
small fraction of the configurational degrees of freedom, \l/ is not the energy (or enthalpy) 
obtained by calorimetric measurements. Therefore, this picture is different from the one 
presented in QQ , where we did compare \1/ qualitatively to calorimetric data. |l8l . Il9l . |20j In 
either case, \1/ must be proportional to the density of STZ's, and must have the dimensions 
of energy per unit volume; therefore it is convenient to write it in the form: 



^(A, A) = p,e Aip(m) 



A 

m = — . 
A 



(2.9) 



The last term on the right-hand side of Eq. ()2.8jl . i.e. Q, is the energy dissipation rate 
per unit volume. The central hypothesis of is that T is simply proportional to the rate 
of energy dissipation per STZ. That is, 



Q(S,A,A) 



To 



AT(s, A,m). 



(2.10) 



We then can use Eqs. fl2.6j) and ()2.7j) to write Eq. (j2.8)l in the form: 



2 C(s) A (T (s) -m)s= U(m) - m if>'(m)) (r(s, A, m) + p(T)) (e" 1/x - a) 

+ if;'(m) A \2C(s) (T(s) - m) - (V(s, A,m) + p(T)) ml + AT(s, A,m), (2.11) 



which can be solved easily for V or, more conveniently, for r = T + p 

r(s, A, m) + p{T) = A ^ j-t 

A -mip'(m) e- 1 /* + ip{m) fe _1 /x - A 

Our equations of motion are now: 



f(5,A,m). (2.12) 



C(S) A (T( 



s) — m 



m ~ 



r m = 2 C{s) (T(s) - mj - ^ T{§, A, m) e~ 1/x ; 



and 



T 



A = f(5,A,m) (e~ 1/x -A): 



(2.13) 
(2.14) 

(2.15) 



Next we must specify the rate factors R(s) and p(T). In |2|, we wrote the dilational rate 
factor in the form 

em^f_^t\ ( , 16) 

r To V v f( T ) J 
where po is a dimensionless prefactor, AVq i1 is the activation volume required for a dilational 
rearrangement, and Vf(T) is usually identified as the free volume. In fact, in we treated 
Vf(T) as a phenomenological function, not necessarily the same as the free volume, and 
evaluated it directly from the measured viscosities with no use of the Vogel-Fulcher or 



Cohen-Grest formulas. 



a, 



2l\ We then suggested, in analogy to the fully nonlinear STZ model 



that the shear rates ought to have the form 

/ A y shear f-x\ \ 

R(s) = exp } } ; AV shear (s) = AV shear e~% (2.17) 

V V A T ) J 

Setting the exponential prefactor equal to unity in Eq. (j2.17j) defines Tq 1 to be the value of the 
dimensional rate R(s)/tq in the limit s — > oo. Since the dominant temperature dependence 
of this rate should occur via the function Vf(T) in the exponent, we may expect that tq is 
at most a slowly varying function of T. With the definition AV shear /v f (T) = a(T), we have 

C(s) = exp [—a cosh(s)] cosh [a sinh(s)] ; (2-18) 

and 

T(s) = tanh [a sinh(s)] . (2.19) 

In the applications to be considered here, there seems to be no reason to expect two different 
activation volumes or two different time constants for dilational and shear rearrangements; 
therefore I shall assume that AVq i1 = AV s ' iear and p = 1. Then 

p(T) = e- a{T) . (2.20) 

Eventually, we shall need to include pressure dependence in the function a(T); but that too 
will be unnecessary for present purposes. 

The final step in deriving equations of motion for A and m is to choose ip(m) so that the 
numerator in the expression for T in Eq. (j2.12j) is non-negative for all values of s. To do this, 
compute the inverse function of T; that is, find the function £ (m) such that = m. Then, 
because T is a monotonically increasing function of its argument, the choice ip'(m) = £(m) 
assures us that both T(s) — m and s — £(m) change sign at the same value of s, and therefore 



that the product of these two factors is never negative. For the specific choice of T given in 
Eq.flTIHD, 



£(to) = In 



A a' 



. 2 /l + m\ 1 , (\ + m 
In ( + — In 



m 



2a 



m 



and 



ip(m) = ip(0) + / £(m) dm. 
Jo 



(2.21) 



(2.22) 



ip(0) is an as-yet undetermined constant which, as it turns out, we shall not need to evaluate. 
The result is 



f (3, A, m) 



A 



M(A, m) 



2C(s) (T(S) -m)(s- £(m)) + p(T) 



where 



M(A, m) = A - m £(m) e~ 1/x + ijj(m) (e~ 1/x - A) . 



(2.23) 



(2.24) 



Positivity of V requires that M(A, m) remain positive along all the system trajectories 
determined by Eqs. tj2.14j) and (j2.15j) in the space of variables A and m. This happens 
automatically so long as all the trajectories start at points where M(A,m) > 0. The locus 
of points along which M(A, m) changes sign is a dynamical boundary for these trajectories; 
the dissipation rate diverges at that boundary, and the trajectories are strongly repelled from 
it in a way that does not allow them to cross into unphysical regions where the dissipation 
rate is negative. An interesting feature of this fully nonlinear case is that the boundary 
always occurs when m is slightly smaller than unity because of the weak divergence of the 
function when m — > 1. For example, setting A = e~ x l x and using Eq. (j2.21j) with a = 2, 
I find that the upper limit of m is 0.983732. We shall see that the interesting values of a 
are generally much larger than this, of order 10 or more, in which case the upper limit of m 
is practically indistinguishable from unity. 

These equations of motion simplify greatly if we note that A = e~ 1//x is always the 
only stable stationary solution of Eq. (j2.15j) . and use this relation to eliminate A from the 
beginning of the analysis. (There seems to be no conventional experimental method for 
adjusting A and x independently of one another.) Then we have 

t e pl = e e- 1 '* C(S) (T(s) - m) ; (2.25) 

_ 2C(s) (T(s) -m) (1-ms) -mp(T) 

T o m — ; 77 — \ ! (2.2o) 

1 — m^(m) 



and 



r 



2C(s) (T(5) - m) (g - £(m)) + p(T) 



(2.27) 



1 — m£(ra) 



Note that the effective temperature \ now appears explicitly only in the strain-rate equation, 
Eq. (j2.25j) . where it determines the density of STZ's that must appear in front of the rate 
factor. Conveniently, the as-yet undetermined energy if)(Q) disappears entirely when we 
assume that A is always in equilibrium with the configurational degrees of freedom at their 
effective temperature. Note also that this approximation will have no effect on any of the 
steady-state calculations presented below. 

Equations ()2.25j) and (|2.2(p l> describe the same exhange of stability at a yield stress that 
we found in earlier papers [l|, |2j, At low temperatures, where p(T) — > 0, the steady-state 
solutions of Eq. (|2.2fi|) are m = T(5) (the jammed state with e pl = 0), and m — 1/5 (the 
flowing state with e pl 7^ 0). The two curves cross at 5 = s y where s y T(s y ) = 1. The jammed 
state is dynamically stable for 5 < s y , and the flowing state is stable for s > s y . For values 
of Oi{T) appreciably larger than unity, the solution of this equation is s y « 1; thus the yield 
stress is s y ~ p,. 

For nonzero p(T), the stable branch of the steady-state solutions of Eq. (|2.2fij) . say m = 

m (s), is 



This function is shown in Fig.(l) along with graphs of m — T{s) and m — 1/5, which are the 
asymptotic values of mo (5) in the limit p — > 0, below and above the yield stress respectively. 
In order that these two sets of curves not lie exactly on top of each other, I have chosen 
a = 6, which will turn out to correspond to a relatively high temperature of about 730 K, 
and have used Eq. (j2.2()j) to evaluate p in Eq. (j2.28j) . 

III. EFFECTIVE TEMPERATURE 

We now need an equation of motion for the effective temperature %■ Before writing such 
an equation, however, it will be useful to discuss some underlying concepts. 





(2.28) 
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The theory as described so far contains three distinct time scales. The first of these is To, 
the roughly temperature-independent time associated with STZ transitions that are driven 
by stresses of order the yield stress or larger. This time will turn out to be very short, of 
order microseconds. A second time scale is r /p(T) = tt, which is the strongly temperature- 
dependent time associated with spontaneous (i.e. stress-independent) thermally activated 
molecular rearrangements. At low temperatures, Tt becomes very long. The third time 
scale is the inverse of the strain rate, (e p ') _1 = r eps , which is determined by the externally 
imposed loading. 

We argued in Q] that r eps is the only relevant time scale for behavior in the regime where 
tq <C r eps <C Tt- This is the situation in which the temperature T is so low that molecular 
rearrangements are not thermally activated, and the actual rearrangements, when they 
occur, are effectively instantaneous. Thus, under steady-state conditions, the number of 
events in which rearrangements occur is not proportional to the time but to the strain. As 
already stated in the Introduction, steady-state deformation with molecular rearrangements 
must produce a steady state of disorder in an amorphous system - a statistical distribution of 
density fluctuations which, in turn, ought to maximize an entropy and therefore be described 
by an effective temperature. Therefore, after very long times, and so long as r eps -C Tt, the 
quantity \ must approach a definite value, say Xoo- A more mathematically precise way of 
saying this, which reminds us that we must be looking in the limit in which both r eps and 
tt are very much longer than To, is that x Xoo if we take the limit e pl — > after T — > 0. 

To make a rough estimate for Xoo, we can use the Stokes-Einstein fluctuation-dissipation 
relation in a way that will require careful discussion. An oversimplified derivation starts by 
noting that, because there is a yield stress s y in these models, the viscosity is r\ = s y / (2 e pl ). 
Then, if the only relevant time scale in the model is T eps , it follows that the diffusion constant 
D, measured over times much longer than r epS) must be proportional to £ 2 e pl , where i is the 
characteristic displacement of a molecule during an STZ-like rearrangement, i.e., roughly 
a molecular spacing. Finally, the Stokes-Einstein relation says that D oc ksT^jrii^ where 
Too = EzXoo/kB- It follows that fesT^ cx £ 3 s y , independent of e pl . 

One problem with this analysis is that the viscosity in the Stokes-Einstein formula is the 
linear response coefficient relating flow to driving force in the limit of vanishing stress and 
strain rate, whereas it is used here at the yield stress. A related and possibly mitigating 
problem is that the "temperature" that we supposedly are evaluating with the Stokes- 
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Einstein formula is a very low-frequency (essentially static) noise strength that determines 
the spatial distribution of energy and density fluctuations but not the rates at which those 
fluctuations vary in time. The temporal rates are determined by the strain rate e pl , which, 
in this limit, is small but much faster than the essentially negligable rearrangement rates 
induced by true thermal fluctuations. 

It helps to visualize the situation as follows. At zero T, for small strain rates, the graph 
of stress as a function of strain-rate consists simply of a horizontal line at s — s y . Now add 
to this system a slow noise source with characteristic frequencies of order, say, uj noise , which 
couples only to the configurational degrees of freedom. Let the strength of this noise source 
be determined by an effective temperature T^. At non-zero T^, the stress versus strain-rate 
graph must start at the origin and rise with a slope 2?7(T 00 ), where ^(T^) is the viscosity 
measured by averaging the stress over times longer than oj~l ise . This section of the curve at 
small i pl will extrapolate to s = s y at a strain rate, say, e^iT^) = s y /(2r](T 00 )). For strain 
rates e pl > ^{T^), the curve returns to s = s y . The Stokes-Einstein relation pertains to the 
portion of this function that goes into the origin at small strain rates. Thus, so long as we 
measure diffusion and viscosity over times longer than co~^ ise , we can estimate the diffusion 
constant: D^T^) oc £ 2 e pZ (Too) and conclude that hsT^ OC £ Sy dbS before, independent of 
what value of e pl = e pZ (Too) we chose at the beginning. 

We can use the above argument - not much more than a dimensional analysis - to make 
an order-of-magnitude estimate for T^. The tensile yield stress that we used in was 
1.9 GPa, which gives us an approximate value for s y . Then ~ 50£ A K, where I a is 
the molecular length scale i measured in Angstroms. If £a ~ 3, then Tqq ~ 10 3 K. A 
similar rough estimate emerges if we guess on dimensional grounds that Ez ~ /i£ 3 , where 
/i is the shear modulus. Then, using the value of Young's modulus given in |^|, we have 
Xoo ~ s y /fi ~ 0.02. If E z ~ 2ev, then again we find ~ 10 3 K. These estimates are 
consistent with Andrea Liu's suggestion [22j that is the glass temperature, which also is 
roughly of the order of 10 3 K for these materials. 

With this understanding of the role and approximate magnitude of the effective tempera- 
ture, we now can deduce an equation of motion for it by returning to the principle of energy 
balance. As noted above, one of the main differences between this model and that discussed 
in Q is that, he,,, the energy ,„ red in the STZ, is only a very small fraction of the energy 
contained in the configurational degrees of freedom. Thus we can assume that the energy 
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dissipated by the STZ's during plastic deformation simply adds to the energy of configura- 
tional disorder. It then seems reasonable to assume that, over the range of temperatures of 
interest here (approximately 550 K - 700 K for the data reported in j^|), the specific heat 
of the configurational degrees of freedom is a constant, say, Cd = ksCo/i 3 , where c is a 
dimensionless number of order unity. The associated configurational energy is Cz>T e ff; and 
the energy-balance equation, i.e. the equation for the rate at which this energy is changing 
per unit time, Cr>T e ^, becomes the equation of motion for T e ff. 

In accord with the discussion in the preceding paragraphs, I propose to write this equation 
of motion in the form: 

C D f eff = Q (l - ^fl) + 6 K( X ) ^k B (T- T eff ). (3.1) 

The first term on the right-hand side says that the energy dissipated during plastic deforma- 
tion, at rate Q per unit volume, is absorbed by the configurational degrees of freedom. The 
second term proportional to Q is the one that says that this process must drive the system 
toward a limiting state of disorder in which T e ff — > T^. This equation can be used only when 
time variations are very much slower than the microscopic rate r _1 , because the preceding 
argument for a limiting value of T e ff is valid only in those circumstances. Thus, we cannot 
expect this theory to be accurate for strain rates higher than about (e /r ) exp (— 1/Xoo)', 
however, such rates are well above the experimental range. 

The term in Eq. (j3.1|) proportional to p{T) says that T e ff — > T in the absence of external 
driving, and does so at a rate which becomes very small at low temperatures. K(x) is an 
equilibration coefficient, defined with a factor eo for convenience. The x dependence of K 
reflects the fact that the equilibration rate must depend on the state of disorder. In what 
follows, I shall assume that 

Kix)=Ke-N*, (3.2) 

so that eo K(x) is proportional to the density of of sites at which the equilibration transitions 
take place. It will be simplest at first to let the equilibration parameter (3 = 1, which means 
that the latter sites are the same as the STZ's. However, as we shall see, other possibilities 
are interesting. 

Next, convert Eq. (j3.1j) into an equation for x by writing, as in Eq. (j2.10|) . 

Q(s, A, A) = ^ e- 1/x r(S, m). (3.3) 

To 
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The function T (not T + p) is 

2C(s)(T(s)-m)(s-£(m)) + p(T)mt(m) 

T(s,m) = * J -± — . 3.4 

1 — m((m) 

Note that the term proportional to p(T) in T disappears in an undriven system because 
m — > in that case. We then find 

^X = e- 1/x r(S,m)( X oo-x) + ^p(T)e-^ ( - X ) ■ (3.5) 

In order to avoid adding another arbitrary constant of order unity, I have used Xoo — fi^ 3 /Ez 
in evaluating the coefficient of T in Eq. ()3.5|) . Equation (|3.5|) . along with Eqs. ()2.25|) . ()2.26j) . 
and ()3.4j) . provides a complete specification of the equations of motion for this model. 



IV. LIMITING BEHAVIORS AT SMALL STRESS 



At this point in the development, it is necessary to rewrite the two-dimensional STZ 
equations of motion in a form in which they can be applied directly to three-dimensional 
experiments, especially those reported in |4|. To do this, I assume that the stresses and 
strain rates are uniform throughout the experimental samples, and follow j^J and ^1 by 
assuming that I can simply replace the variables s, e pl , and m by traceless symmetric tensors. 
In the case of a uniform sample with uniaxial applied stress in, say, the x direction, and 
free, stressless surfaces normal to the y and z axes, each of these tensors is diagonal with 
elements proportional to (1, —1/2, —1/2). The total stress tensor has only one nonzero 
element, a xx = a. Define m 2 = (1/2) m^- m^- = (3/4) m 2 x , so that m = ^j3/4m xx ; and, 
similarly, s 2 = (1/2) % % = (3/4) s 2 xx1 s = ^3/4 s xx . The only way in which this analysis 
differs from that in j2| is that, here, we must form tensorial generalizations of the functions 
T(s) and £(m). This can be done most simply by writing %j(s) = (%/|s|) T(s), and 
£ y (m) = (mij/lmD^ijim). 

With these transformations, we recover precisely our earlier formulas, Eqs. ()2.25J) . (|2.26j) . 
f!3.4|) . and ()3.5|) as the xx components of the tensor equations. The single difference is that, 
because the experimental strain rate e p x l x is not rescaled as are s and m, the parameter eo 
in Eq. ()2.25j) is replaced by e' = ^4/3 eo- The low-temperature exchange of stability still 
occurs when s = s y where s y T(s y ) = 1; therefore, at the temperatures of interest here, we 
still have s y = 1, s y = p. The experimental data is expressed in terms of the tensile stress, 
which becomes a = (3/2) s xx = \f2>ps = a y s, where a y is the tensile yield stress. 
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The first quantity that we must compute is the Newtonian viscosity t/jv, that is, the linear 
viscosity in the limit of vanishingly small stress and strain rate. As in j^J, comparing our the- 
oretical 777V with the experimental measurements reported in |4] provides initial constraints 
on several of the parameters that appear in our equations. In the small- s limit, we have 
T{s) ~ as, and C(s) ~ C(0) = exp(— a). Then we can deduce from Eqs. (|2.19|) . (|2.2U|) and 
f!2.28|) that, to lowest (linear) order in s, 

2 _ 



m (s) 



a s. 



Using Eq. ()2.26j) . and noting that the product mo(s) s is small of order s 2 , we have 



2C(s) (T(s)-mo(C 
so that, using Eq. (j2.25|) . we find 



-e- 1/x C{s) (T(s) 

T~0 V 



mo(s)p(T); 



2 r 



m (s) p(T). 



(4.1) 



(4.2) 



(4.3) 



To evaluate x from Eq. (j3.5j) . note from Eq. ()3.4|l that T is small of order s 2 , so that \ ~ T/Tz 
in the small-s, steady-state limit. Therefore, 

e' 



— a(T) s exp 
3r 



' T 



and 



r) N {T) = lim 



V3pr 



cxp 



a(T) 



T 7 



(4.4) 



(4.5) 



,0 2eS4 ^a(T) 

We also can use this analysis to compute the temperature-dependent stress relaxation 
rates discussed in j^. In these measurements, samples first were compressed at relatively 
small strain rates and then held at fixed total strain e total while the stress was measured as 
a function of time. The equation of motion that we must solve therefore is 



•total 

XX 



0. 



(4.6) 



where E is Young's modulus. Using the preceding small-stress approximations, and again 
assuming that x is thermalized in these experiments (x ~ T/Tz), we find 

s 



e' E 

■ ct(T) exp 



T 7 



tr(TY 



which implies that the exponential relaxation time is 

Vn(T) 



E 



(4.7) 



(4.8) 
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This relation is consistent with the conclusion of Lu et al. [4] that both t]n{T) and t r [T) 



■ Q 



sea 
by 



e with the same temperature dependent rate factor. However, if E = 96 GPa as reported 



,1 



these theoretical values of t r are too small by a factor of about 50. This discrepancy 
may be due to the thermalization assumption, which might not be consistent with the way 
in which these measurements were made. 

In we used our expression for r?jv(T) and the experimental values for this quantity 
given in j^J to obtain estimates of p(T) (up to a scale factor) for eight separate values of 
the temperature. However, our present formula for the Newtonian viscosity, Eq. (j4.5|) . is 
more complicated than the one in because it now contains a temperature-dependent STZ 
density, proportional to exp (— 1/x) ~ ex P (~Tz/T), as well as the temperature-dependent 
rate factor p(T). Moreover, it will be useful for present purposes to have a smooth func- 
tional representation of p(T) rather than just values at separate points. Accordingly, I 



have fit Eq. (j4.5|) to the experimental data using the Cohen-Grest formula 21] as a purely 



phenomenological fitting function for a(T). That is, in Eq. ()2.20j) . I have used 

a(T) = Tr (4.9) 

T-To + JiT-ny + TtT 



where Tr, To, and T\ are fitting parameters with the dimensions of temperature. 

Clearly, we cannot obtain a unique fit for all the parameters in Eqs. (j4.5)l and ()4.9)1 from 
just the viscosity data, so we now must make some physical assumptions. The guiding 
principle is to make the simplest possible choices, and to add complications only if they 
become necessary. In this spirit, we may assume that there is only one temperature that 
characterizes the glass transition. In Eq. ()4.9|) . that temperature is Tq. If we then adopt 
Liu's hypothesis [22j that T^ is the glass temperature, we should choose T^ = T . On the 
basis of various clues, including calorimetric analyses, I estimate that T^ = T = 800 K. 



This value is consistent with my guess that the temperatures used in the experiments of 
are all well below T ; that is, the behaviors seen in these experiments seem characteristic of 
states in which the material is softening rapidly with increasing T but is still stiff enough to 
resemble a solid in resisting deformation. 

Next we must estimate Tz- In the preceding dimensional analysis, we guessed that 
Tz/Too ~ /■*/% ~ 50, which would imply that Tz ~ 40, 000 K. This estimate, however, is 
uncertain by a least a factor 2. A better strategy, I think, is to assume that the Newto- 
nian viscosity is dominated at the higher temperatures shown in Fig. 10, by the factor 
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exp (Tz/T) in Eq. (j4.5|) . That fitting strategy yields Tz = 25, 000 K, which is within the 
previous uncertainty and is the value that I will use here. It means that Ez — 2 ev, which 
seems plausible for a vacancy-formation energy. 

Figure (2) shows the fit to the Newtonian viscosity as a function of temperature with 
T z = 25, 000 if, T = 800 K, T R = 600 K, T x = 28 K, and rj = fiT /e = 2 x 10- n Pasec. 
The eight points at the lowest temperatures in Fig. (2) are the ones that we used in Q; the 
four points at higher temperatures are also taken from j^J, Fig. (10). Figure (3) shows the 
corresponding function a(T). The values of a in the range of experimental interest, roughly 
600 K to 700 K, are of order 8 to 15, that is, in about the same range as the values that 

n 

Falk and I M found to fit the original MD simuations. 

n 

As in |2|, I assume that the tensile yield stress at the experimental temperatures is the 
same as the room temperature value reported in j^J, i.e a y = 1.9 GPa. Thus fta y /\^3 = 1.1 
GPa. With the above value of r] Q , we have e /ro = 6.3 x 10 19 sec _1 . 

V. ANALYSIS AND COMPARISON WITH EXPERIMENTS 

We are now ready to explore the properties and experimental predictions of this effective 
temperature theory at values of the stress and strain rate where the response to loading 
becomes nonlinear. Look first at the steady-state solutions obtained by using Eq. (|2.25|) 
to compute the strain rate (with €q — > e ), and by setting the time derivatives on the left- 
hand sides of Eqs. (|2.26|) and ()3.5|) to zero. The stresses and strain rates found in this 



way correspond to those obtained by Lu et al. 4] from the late, steady-state stages of 
their constant-strain-rate measurements. The steady-state values of the reduced effective 
temperature \ may, in principle, be obtained by calorimeteric measurements as discussed 
below. 

It is simplest to start by setting (3 — 1 in Eq. (j3.2)l which, as mentioned earlier, means 
that the thermal fluctuations that drive the effective disorder temperature toward the tem- 
perature of the heat bath occur predominantly at the STZ sites. The only other adjustable 
parameter in steady state is k. Figure (4) shows the dimensionless stress s as a function 
of the scaled strain rate r/jv(T) e^ x for /3 = 1, k — 2, and for four different temperatures 
T in the range of the metallic glass data |4j. It should be compared with Fig. (3) of [2j, in 
which these curves lie accurately on top of one another up through the yield stress. This 
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scaling behavior, which was discovered experimentally by Kato et al. 3] and explored in 
more detail by Lu et al. 4|, clearly is broken here. The trend toward lower stresses at lower 
temperatures can be understood as a nonlinear property of the effective temperature theory. 
As the strain rate increases, T e ff increases, and the driving force required to maintain that 
strain rate decreases accordingly. Because 77jv(T) increases rapidly with decreasing T, we 
may understand the curves that are plotted as functions of tjn{T) e^ l x in Fig. (4) to be a 
sequence in which T e ff increases as T decreases. 

The important question is whether this non-scaling behavior is ruled out by experiment. 
Figure (5) shows a direct comparison of the data from [4] with theoretical curves for eight 
different temperatures, in analogy to Fig. (5a) in Q]. The value k = 2 was chosen to optimize 
the fit to the data at 643 K. If we use the temperatures cited in |4|, the agreement is 
reasonably accurate for 623 K and above (apart from a few apparently outlying points), 
and also (perhaps fortuitously) is satisfactory for the two points at 573 K. The theoretical 
curves for the latter set of temperatures are shown by solid lines in the figure. 

However, the agreement is not so good at 593 K and 603 K , and is especially poor at 
613 .ft". In interpreting this disagreement, remember that we evaluated p(T) in |2j point by 
point from the viscosity data given in J] , and then checked that these values were the same 
as those that the latter authors had used in scaling their strain rates. Thus we did not use 
the nominal values of the temperature T in any of those analyses. Here, on the other hand, I 
have fit p(T) by an analytic expression, Eqs. (j2.20j) and (|4.9j) . and have used this function of T 
in plotting the curves shown in Fig. (5). The problem is that the viscosity data for the lowest 
four temperatures shown in Fig. (2) does not fit onto a smooth curve. Since there is no reason 
;o believe that the material is undergoing any qualitative change in this temperature range 

. we must presume that either the reported temperatures or the Newtonian viscosities 
- or both - are inaccurate. Accordingly, instead of using the nominal temperatures 593 K, 
603 -ft", and 613 K in drawing the theoretical curves in Fig. (5), I have used 588 -ft", 598 -ft", 
and 607 -ft" respectively, and have indicated these results by dashed lines. Note that these 
small shifts in temperature produce large shifts in the predicted Newtonian viscosities, and 
thus move the data points closer to the smooth curve in Fig. (2). More importantly, the low- 
temperature data is largely in the region where the linear Newtonian behavior is becoming 
nonlinear superplasticity; therefore the ability of the theory to account for the data seems 
significant. In short, although the effective temperature theory systematically departs from 
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the strong scaling behavior obtained in [2j], it appears that these departures are within the 
present uncertainties in the experimental data. 

Our one remaining point of contact with the data of Lu et al. ^| is the transient response 
shown in their constant strain-rate experiments. As in we can use these experimental 
results to obtain separate estimates of the parameters eo and To instead of just their ratio. 
To compute the corresponding stress-strain curves, write the equation of motion for the total 
strain e total (including elastic strain) in the form 



and solve this simultaneously with Eqs. ()2.26|) and (|3.5J1 for s, m, and \ at fixed e 1 "^ 1 . In 
preparation for plotting stress-strain curves, we can let e^ aZ replace time as the independent 
variable, in which case e' appears separately as well as in the combination co/Vo Thus, fitting 
the transient response yields separate estimates for e' Q and r . 

The parameter e' must be a large number in this version of the STZ theory, because 
the fraction of the volume covered by STZ's is proportional to e' Q exp (— 1/x), not just to 
e' by itself as in |2(. At x — Xoo, this fraction would be unity if e' were about 3 x 10 14 . I 
find that e' = 10 14 works well for making the theoretical stress-strain curves agree with the 
experimental ones shown in |4J], Figs.(l) and (2). With that value, the equilibrated fractional 
density, e' exp (—Tz/T) is of order .002 at T = 648 K; thus the effective temperature theory 
produces estimates of the STZ density that are in accord with the idea that this density 
should be small, which was not necessarily the case in Q]. With e' = 10 14 , and our earlier 
estimate e' Q /r = 6.3 x 10 19 sec _1 , we have r ~ 10~ 6 sec, which seems reasonable if we 
remember that Tq 1 is the STZ transformation rate in the limit of infinite applied stress. 

With these parameters, plus k = 2, (5 = 1 and Co = 1 in Eq. (j3.5j) . and E/a y = 50, the 
stress-strain curves are essentially identical to those shown in Figs. (1) and (2) of which, 
in turn, were similar to the experimental ones shown in P|. In this fully nonlinear theory, 
the initial rise of the stress is accurately determined by Young's modulus, instead of being 
too small because the plastic response was unrealistically enhanced at small stresses in the 
quasilinear version. Typical stress-strain curves, analogous to those in j2|, are shown in 
Figs. (6) and (7). As in p, initial values of s and m are zero. I have assumed that the initial 
value of x is always equal to T/Tz, that is, that the samples are completely equilibrated 
initially by annealing at the experimental temperatures. 




(5.1) 
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Turn now to the thermodynamic properties of this theory, which were missing in the 
earlier version jsj but now can be explored in detail. Given our assumption that the con- 
figurational energy is simply proportional to the effective temperature, we can convert the 
equation of motion for \ in the absence of driving, i.e. Eq. (J3.5j) with T = 0, into an equation 
for the specific heat measured in a DSC experiment. Let the heating rate be h = T. Then 
(J3.5|) becomes 

* = jL«£2l e Hv x fT_ y (5 . 2) 

dT c t h \T Z J 

The left-hand side is equal to the specific heat in units CdT z . Figure (8) shows dx/dT 

computed by solving Eq. (j5.2j) with k = 2, (3 = 1 and Cq = 1, and a heating rate of 10 K per 

minute. The initial temperature used for integrating this equation was 550 K\ but the results 

are insensitive to this value so long as it is low enough. The different states of the system are 

specified by the initial values of x which, in temperature units (i.e. expressed as T e ff = Tz x) 

were chosen here to be 630 K, 640 K , 650 K , and 660 K . These curves resemble those found, 

ta exa m ple, In De Hey et a.. Q „ r M „ Str a et a.. Q Eventually they S nou ld be eo m pa r e d 

with data for Zr 41 2 Ti 13 8 CU12.5 Ni 10 Be 2 2.5 such as that shown in Fig. (2) of Busch et al. 

but that analysis would best be carried out in connection with experiments like those 

discussed in the next paragraph. Not shown in Fig. (8) is the prediction from Eq. (J5.2j) that, 

when the system is fully annealed at lower temperatures so that the initial T e jf is well below 

600 K, the spike becomes very sharp and moves to temperatures above 700 K. 

The differences between the areas under specific heat curves of the kind shown in Fig. (8) 

are equal to the differences between the enthalpies of systems with the corresponding initial 

values of x- Those values can be controlled experimentally by shearing the system at fixed 

bath temperatures and strain rates for long enough times that they achieve steady state. 

On the theoretical side, we can compute the values of x as functions of temperature and 

strain rate by finding the steady-state solutions of Eqs. ()2.26|) and ()3.5|) . as we have done 

to obtain the steady-state stresses in Figs. (4) and (5). Thus our steady-state values of x 

can be determined experimentally. Precisely such measurements have been performed by 

De Hey et al. who interpreted their function x as a reduced free volume instead of 

a reduced effective temperature. The two interpretations may be effectively equivalent for 

systems held at constant pressure because, under that condition, the change in volume will 



be proportional to the effective temperature. 

Figure (4) in De Hey et al. [3| shows values of x as functions of strain rate at three 
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different temperatures for thin ribbons of amorphous PcUo Ni 40 P 2 o- Without detailed in- 
formation about other parameters of the kind obtained here from J4J], we cannot try to 
reproduce the results of |3| theoretically. Instead, I have used the parameters determined 
here for amorphous bulk Zr 4L2 Ti 138 CU12.5 Ni 10 Be 2 2.5 to compute a graph analogous to the 
one in [3]. The results are shown in Fig. (9) for four different bath temperatures T. Note 
that x approaches T/Tz at low strain rates, and goes to Xoo = 0.032 at large strain rates. 



At intermediate rates, such as those shown in [18j], the values of x decrease as T increases, 
consistent with the idea that the number of STZ's needed to sustain a fixed strain rate de- 
creases when thermal fluctuations assist the transitions. As anticipated in 3|, these curves 
cross each other as they move to small strain rates. Figure (9) implies that this crossover 
might be observed experimentally in Zr 412 Ti 13 8 CU12.5 Ni 10 Be 22 . 5 . 

All of the preceding calculations have been based on the equation of motion for x, 
Eq. ()3.5p . with the equilibration parameter f3 set equal to unity. Remember that the quantity 
/SksTz is a characteristic formation energy for configurational fluctuations that drive the 
effective temperature T e jf toward the bath temperature T. A value of f3 smaller than unity 
implies that these fluctuations occur more frequently than the STZ's, which seems plausi- 
ble. (The opposite situation, (3 > 1, might also occur.) Figure (10) shows what happens 
to the steady-state stress versus strain rate curves if we choose (3 — 0.5. In order to be at 
least roughly consistent with experimental data, that is, in order that the two terms on the 
right-hand side of Eq. ()3.5j) be of comparable size when x is near Xoo, we must choose a much 
smaller value of k than previously. Specifically, k = 10 -8 for the graphs shown in Fig. (10). 
The most important new feature is that the curve for the lowest of the four temperatures 
shown here, T = 573 K no longer remains below the others as it does in Fig. (4), but now 
rises above and goes through a maximum and then a minimum before returning to approxi- 
mately its previous behavior. This multi-valued property is seen more clearly if, instead, we 
plot the reduced effective temperature x as a function of the stress, as shown in Fig. (11). 



The multi-valued behavior o: 
analyses by Olmsted et al. [26 



Y implies a shear-banding instability. (See, for example, the 
22|,|2S| of shear banding in several similar situations.) In the 



usual simple-shear experiment in a strip geometry, the shear stress remains constant across 
the sample in order to satisfy force-balance. If the externally imposed shear rate is chosen 
so that the stress lies in the multi- valued region, then the sample will have to break up into 
regions of large and small flow or, equivalently, high and low effective temperature. That 
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is, the system will encounter a shear-banding instability and most likely will fail via shear 
fracture. Figure (11) indicates that this instability appears only at temperatures lower than 
about 648 K, and that the onset stress increases with decreasing temperature. The figure 
also implies that, even at low temperatures, uniform flow should be stable at sufficiently 
large driving forces. 

A satisfactory theory of shear banding also needs a length scale, because it must de- 
scribe a smooth transition between the jammed and flowing regions of the material. The 
effective temperature theory suggests that a natural way to introduce this length is to add 
a diffusion term proportional to V 2 x t° t ne right hand side of Eq. (|H.5j) . The associated 
diffusion constant will be very much smaller than the ordinary thermal diffusion constant 
because configurational disorder must diffuse extremely slowly at temperatures below the 
glass transition. We may even be able to estimate the magnitude of this diffusion constant 
from the arguments presented in Sec. III. Thus, the effective temperature theory seems to 
be giving us a clue about how to solve the long-standing problem of identifying an intrinsic 
length scale for shear localization. A fully detailed development of these ideas, however, is 
beyond the scope of the present paper. 

VI. CONCLUDING REMARKS 

My main conclusion is that the effective-temperature version of STZ theory looks promis- 
ing but is far from being quantitatively confirmed by comparison with experimental data. I 
see several directions for future investigations. 

First, there is a need to combine mechanical and calorimetric measurements, in the 
manner described by De Hey et al. Q|, in order to test predictions of the kind shown 
in Fig. (9). Such experiments may come as close as is possible to actually measuring the 
effective temperature and learning whether it behaves as predicted. It also would be useful 
to repeat the mechanical experiments with enough precision to test the predicted deviations 
from scaling shown in Figs. (4) or (10). For this purpose, it might be well to use other 
materials such as polymeric glasses or, perhaps, colloidal suspensions in order to control the 
experimental conditions more precisely than seems possible with amorphous metals. 

A second direction for further research is to develop the theory of shear banding along 
the lines described above and to test the results experimentally. For example, it should be 
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possible to predict and measure the onset of spatial instability as a function of temperature 
and applied stress. 

In my opinion, the principal theoretical question left unanswered is the temperature 
dependence of the rate factor p(T). This factor has been determined empirically here from 
experimental measurements of the Newtonian viscosity, with no theoretical justification 
whatsoever. The STZ theory described in this paper and in |2j differs most markedly from 
the earlier flow-defect theories in that we ascribe the non-Arrhenius behavior of the viscosity, 
not to the density of defects, but to the transition rates. Thus the kind of analysis that 
was used in deriving the Cohen-Grest formula [2]] seems unlikely to be applicable. Instead, 
we must learn how to perform a truly nonequilibrium analysis of the processes by which 
configurational rearrangements occur, both spontaneously and as driven by imposed stresses. 
This is a large challenge - perhaps the same as that of understanding the fundamental nature 
of the glass transition itself. 
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FIG. 1: Graph of the function itiq(s), s = s/ fi ~ s / s yi f° r a = 6. Also shown are the related 
functions T{s) and 1/s. 
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FIG. 2: Experimental values of the Newtonian viscosity t]n(T) taken from Lu et al.,|J] and the 
analytic fit to these points obtained by choosing parameters in Eas. H2.2U|) and ()4.9j) as explained 
in the text. 
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FIG. 3: The function a(T) determined by fitting the viscosity data shown in Fig. (2). 
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FIG. 4: Graphs of the scaled deviatoric stress s = s/s y as functions of the scaled, steady-state 
strain rate t]n{T) for four different temperatures. 
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FIG. 5: Tensile stress as a function of the steady-state strain rate e^ x (in sec -1 ). The data points, 
taken from Lu et al. J], correspond to the nominal temperatures as shown, and the solid theoretical 
curves are computed using those temperatures. The dashed curves are computed using, from left 
to right, T = 588 K, 598 K, and 607 K. 
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FIG. 6: Theoretical stress-strain curves for I™™ 1 = 0.1 sec at three different temperatures. 
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FIG. 7: Theoretical stress-strain curves for T = 643 K at four different strain rates e! 




FIG. 8: Scaled specific heat curves, dx/dT, corresponding to simulated DSC measurements at 
a heating rate of 10 K per minute. The initial states, in order of decreasing peak height, have 
effective disorder temperatures T e ff = 630 K, 640 K, 650 K, and 660 K. 
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FIG. 9: Steady-state values of the scaled effective temperature x as functions of strain rate at four 
different temperatures. 
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FIG. 10: Graphs of the scaled deviatoric stress s = s/s y as functions of the scaled strain rate 
t]n(T) e^ x for four different temperatures. This figure is analogous to Fig. (4), but the equilibration 
parameter (5 has been set equal to 0.5. 
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FIG. 11: Graphs of the scaled effective temperature x as functions of the scaled deviatoric stress 
5 = s/sy at four different temperatures. As in Fig.(10), the equilibration parameter is (3 = 0.5. 



29 



